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Abstract 

We report the results of a search for the decay e'^e~ , /u^/u^ and e^/i^ based on an analysis 

of 78 fb~^ of data collected by the Belle detector at KEKB. No candidate events have been found. 
Upper limits on the branching fractions are calculated at the 90% confidence level: B{B^ e+ 
e-) < 1.9 X 10-^, B{B^ n+ n') < 1.6 x lO"'^, and B{B^ e± /x^) < 1.7 x 10"^. A limit on 
the Pati-Salam leptoquark mass Mlq > 46 TeV/c^ is obtained at the 90% confidence level. 

PACS numbers: 13.20.He,14.40.Nd 
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The B meson decays e+e and B^ jji^jj, can proceed at a low rate through the 

Flavor- Changing Neutral Current (FCNC) process. The Standard Model (SM) predictions 
for the branching fractions are (2.34 ±0.33) x lO'^^ and (1.00 ±0.14) x 10"^° for S° e+e" 
and 5° decays, respectively [l|. These branching fractions could be two orders of 

magnitude larger in models with two Higgs doublets and Z- mediated FCNC jH . The decays 
B^ — > are essentially forbidden in the SM by lepton number conservation, apart from 

a vanishingly small probability permitted by neutrino oscillation. They can occur, however, 
in the Pati-Salam Model [3] or supersymmetric (SUSY) models The charge conjugate 
modes are implicitly included throughout this letter. 

In this letter we report the results of a search for the decays 5° —>■ e^e~ , fi~^fj^~ and e^/x^ 
(collectively denoted as 5° — > using 78 fb~^ of data, corresponding to 85 million BB 

pairs, collected with the Belle detector at KEKB , an asymmetric e~^e~ collider operating 
at the T(4S') resonance. The Belle detector is a large-solid-angle magnetic spectrometer 
that consists of a three-layer silicon vertex detector (SVD), a 50- layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters (ACC), time-of-flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron fiux-return located outside of 
the coil is instrumented to detect mesons and to identify muons (KLM). The detector 
is described in detail elsewhere 0. 

For each primary charged track, the impact parameter relative to the run-by-run inter- 
action point is required to be within 2 cm along the z axis (aligned with the positron beam) 
and within 0.08 cm in the plane transverse to this axis, and the transverse momentum must 
exceed 100 MeV/c. At least five charged tracks including non-primary tracks are required. 
For each event the visible energy in the detector is required to be greater than 0.6 GeV. 
Monte Carlo (MC) studies show that the radiative Bhabha events (100%), as well as two- 
photon (94%) and t~^t~ (85%) events, are rejected by these criteria. The two-photon MC 
studies use the generator AAFH 0, which includes all tree-level diagrams to fourth order 
in aQED- 

Electron tracks are distinguished from hadron tracks by examining the amount of energy 
and the transverse shower profile of the associated ECL cluster as well as light yield in 
the ACC. The dE/dx measurements in the CDC is also taken into account. Muon tracks 
are distinguished from hadron tracks by their range and transverse scatteringin the KLM. 
We select electrons and muons by requiring that the normalized likelihood ]a |^ that the 
track is a lepton rather than a hadron exceeds 0.9. The electron and the muon detection 
efficiencies are 88% and 89% which are studied via J/ip ^ e~^e~ and J/ip ^ in the 

same track momentum range, 1.5 GeV/c - 4.5 GeV/c. The corresponding misidentification 
rates from kaon tracks and pion tracks are 1.7% {K /i), 1.0% (vr fx), < 0.005% 
{K — i> e), and < 0.1% (tt e) which are studied using kaons and pions from the decay 
chain D* — >■ D^tt^ — > {K^tt^)tt~ in the same track momentum range, 1.5 GeV/c - 4.5 
GeV/c. 

B meson candidates are formed from pairs of oppositely charged leptons, not necessarily 
of the same flavor, and selected using two kinematic variables deflned in the T(45') center-of- 

mass (CM) frame: the beam-energy constrained mass, Mbc = i/-Ebeam ~ Pb^ energy 
difference, AE = Eb — Ey^eam, where pb and Eb are the momentum and energy, respectively, 
of the B candidate and -Ebeam is the beam energy. The signal region is deflned as 5.27 GeV/c^ 
< Mbc < 5.29 GeV/c^ and \AE\ < 0.05 GeV for all channels. The sideband region, used for 
the background determination, is deflned as 5.20 GeV/c^ < Mbc < 5.26 GeV/c^ and \ AE\ 
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< 0.30 GeV. 



TABLE I: The number of expected events are normalized to the integrated luminosity of the data 
sample. The results show that the backgrounds in MC, off-resonance data and on-resonance data 
are in good agreement after the N > 5 cut. 



Mode Artrk=3,4 ^trk>5 N'^^ll^ 

T+T- 14 ±4 2.4 ±1.7 

two photon 250 ± 109 16 ± 5 

qq{q = u,d,s) 6±1 82 ± 7 1.3 ±0.9 

qq {q = c) 54 ± 6 396 ± 16 3.8 ± 1.5 

off-resonance data 396 ± 46 431 ±46 10 ± 7 

data 441 473 12 

The decays K^tt^ and vr+vr^ do not contribute measurably to the back- 

ground because of their small branching fractions, the small lepton fake rates, and, for the 
K^TT~ mode, the shift of the B candidate out of the AE signal region due to the use of lepton 
masses for the tracks. Background from other B decays is also found to be negligible. Back- 
grounds from continuum e+e~ ^ qq {q = u, d, s, c) production are suppressed using seven 
variables that characterize the event topology: five modified Fox- Wolfram moments jlol |. 
S± and the cosine of the polar angle 6'thr of the B candidate's thrust axis Events 



e~^e~ cc with two semileptonic charm decays produce B^ candidates which are 

not suppressed by lepton identification cuts; without special treatment, these events would 
be the dominant source of background. We use four additional variables to suppress these 
events: the missing energy -Emiss and momentum Pmiss, the cosine of the angle ^pi between the 
missing momentum and the B candidate's thrust axis, the cosine of the angle 6'pb between 
the missing momentum and the thrust axis of the remaining particles in the event. All 



eleven variables are combined into a Fisher discriminant |12[ . We obtain probability density 
functions (PDFs) for the signal and for continuum background as functions of this Fisher 
discriminant (F) and of the cosine of the polar angle (cosOb) of the B candidate's flight 
direction. We then obtain signal and background likelihoods, Cs and Cb, from the product 
of the PDFs for F and cos^^, and demand that CTZ = Cs/{Cs + Cb) exceeds 0.9. This 
criterion was chosen to optimize the performance of the signal significance {Ns / a/ Ns + Nb) 
using both MC and data samples, where Ns and Nb are predicted signal yields and pre- 
dicted background yields, respectively. Assuming the branching fractions of B^ £^i~ are 
10"'' ~ 10~^^, we then obtain the signal yields from the product of the total number of BB 
events, the signal efficiency from MC and the assumed branching fractions. The predicted 
background yields are obtained by fitting the distribution of AE from the sideband data 
projection in the signal region. The backgrounds in the on-resonance sideband region are 
compared to MC and to data taken 50-60 MeV below the T{4S) resonance. The results are 
summarized in Table HI Distributions of -Emiss, cos6b, F and CTZ are shown in Fig. ^ 

The signal efficiency for each mode, after application of all selection criteria, is determined 
by MC simulation and itemized in Table HTl The systematic error of this efficiency arises from 
tracking efficiency (1% per track), electron and muon identification efficiencies (4% and 2% 
per lepton, respectively, determined from a study of J/ip decays), the charged-track-count 
criterion (10%, measured from the uncertainty in the efficiency using a. B^ K~^7t~ control 
sample), and the likelihood ratio cut (4%, based on a study of B~ D^-n^ {K^-n^)-n^). 
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FIG. 1: Open histograms are from signal MC and the points with error bars are from sideband 
data. The arrow in the (d) CTZ plot indicates the optimized selection. 

The total systematic uncertainty is obtained by first summing the correlated errors linearly 
and then adding the non-correlated errors in quadrature. 

Since no signal events are found for any mode, we determine 90% confidence level upper 
limits on the branching fractions by an extension of the Feldman- Cousins method, which 
takes systematic uncertainties into account The qq events are the main background 

after all the event selections. In order to estimate the number of expected background 
events in the signal box, a qq MC sample (corresponding to an integrated luminosity 1.6 
times greater than that of the data) is used to fit the shape of the background distribution 
in AE and Mbc in the sideband region. These distributions are scaled to the number of 
observed data events in the sideband region; the number of expected background events is 
then obtained from the signal-to-sideband ratio. The dominant systematic uncertainty on 
the background estimation comes from the statistical error of the number of sideband data 
events. The results are listed in Table |n] 

TABLE II: Summary of the £~^i~ search, where e is the reconstruction efficiency, A'obs is 

the measured number of events in the signal region, A'^p is the expected background in the signal 
region and BF is the 90% confidence level upper limit for the branching fractions. 



Mode 




e[%] 


A-obs 


^ 'cxp 


BF(10-'^) 






14.3 ±2.0 





0.30 ±0.12 


< 1.9 






16.9 ±2.0 





0.19 ±0.10 


< 1.6 






15.8 ± 1.9 





0.22 ±0.10 


< 1.7 



These upper limits can be used to constrain certain extensions of the Standard Model. 
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For example, the Pati-Salam model ^3] predicts the existence of lepto quarks — heavy spin-1 
gauge bosons that carry both color and lepton quantum number — that mediate the decays 
e^jji^ 3 3- The branching ratio is related to the leptoquark mass Mlq by 



where 



We follow the method in Refs. |1(tI.I17| and use Fb = 175 ± 30 MeV for the decay con- 
stant, niBO = 5279.3 ± 0.7 MeV/cHlJfor the B° meson mass, nib = 4.3 ± 0.2 GeV/c^ 
for the 6-quark mass and rrit = 176.0 ± 6.5 GeV/c^ '03] for the t-quark mass. We use 
(Mz) = 0.117 Hi which is evolved to Mlq using the Marciano approximation jl^, assuming 



no other colored particles exist between and Mlq. We obtain a limit of Mlq > 46 TeV/c^ 
on the mass of the Pati-Salam leptoquark at the 90% confidence level. This bound holds 
in the case where the r lepton is associated with the strange and charm quarks within 
the Pati-Salam theory. Previous lower limits for this case from B^ i^u jisl . l2o| and 
B^ — » /ie 0, 0| are weaker than the bound presented here; the limit from lepton univer- 
sality in TT"*" ^ f+z/ is stronger 3 3- 

In conclusion, we find no evidence for B^ — > e+e^, fJ^^fJ^^ , and e^/i^ decays in 85 million 
BB events, and set upper limits on the branching fractions of 1.9 x 10~^, 1.6 x 10~^, and 
1.7 X 10"^, respectively, at the 90% confidence level. We use the latter upper limit to set a 
lower bound on the mass of the Pati-Salam leptoquark of 46 TeV/c^. 
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